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Abstract. In this article, we investigate the stochastic resonance (SR) effect in a 
finite array of noisy bistable systems with nearest-neighbor coupling driven by a 
weak time-periodic driving force. The array is characterized by a collective vari- 
able. By means of numerical simulations, the signal-to-noise ratio (SNR) and the 
gain are estimated as functions of the noise and the interaction coupling strength. 
A strong enhancement of the SR phenomenon for this collective variable in com- 
parison with SR in single unit bistable systems is observed. Gains larger than 
unity are obtained for some parameter values and multi- frequency driving forces, 
indicating that the system is operating in a non-linear regime albeit the smallness 
of the driving amplitude. The large SNR values observed are basically due to the 
fact that the output fluctuations are small and short lived, in comparison with 
their typical values in a linear regime. A non-monotonic behavior of the SNR with 
the coupling strength is also obtained. 

1 Introduction 

The study of the response of complex systems formed by coupled nonlinear noisy elements sub- 
ject to the action of a time-periodic force brings into the picture a new ingredient with respect 
to the response of a single unit: the effect of the interactions. The interplay of noise, intrinsic 
nonlinear dynamics, driving forces and interactions render the response a very rich stochastic 
process. Certain aspects of the response are associated with what is called stochastic resonance 
(SR) [HE], a phenomenon of interest to scientific fields other than Physics (as discussed, for 
instance, in [3]). SR is characterized by the non-monotonic behavior of some quantifiers with 
the noise strength. An enhancement of SR effects with respect to those observed in single units 
have been reported in several papers [HEHTIIH] . More recently [QjJlOtfTT] , we have analyzed the 
collective response of a finite set of globally coupled bistable systems. We demonstrated that 
the response of the collective variable shows SR effects which are indeed much enhanced with 
respect to those in single bistable units. Actually, in those arrays, SR gains larger than unity 
were observed, we believe for the first time, for subthreshold sinusoidal driving forces. Those 
findings indicated that the arrays were indeed operating in nonlinear regimes. Similar results 
have been reported for single bistable units driven by a suprathreshold sinusoidal force |12j , by 
subthreshold multi-frequency forces [T3lfl4. 15, 16, T71fT8] . or by a subthreshold sinusoidal force 
in the presence of a strong, high-frequency monochromatic signal [19] . Nevertheless, to the best 
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of our knowledge, gains larger than unity have never been observed when the single bistable 
unit is driven by just a subthreshold sinusoidal force. 

In a later work [20] , the response of the collective variable of a finite array of globally coupled 
units to a rather weak driving force was analyzed. Two aspects of the response were studied: 
SR and noise induced phase frequency synchronization. The fact that SR effects were greatly 
amplified and that very good phase frequency synchronization was found in the response to a 
weak driver led us to indicate that the arrays were indeed operating in nonlinear regimes. 

In this work, we extend the results of the previous work to other topologies of the array. 
In [8], arrays of nonlinear bistable units with nearest neighbor coupling are considered. In that 
work the response of a single oscillator is found to be enhanced due to the coupling to the other 
units in the chain. In [21] , the authors consider an array of driven spins with Glauber dynamics 
and nearest neighbor interaction. They study not only the response of a single spin, but also the 
behavior of the global magnetization of the array. Their analysis of the system response is based 
on the linear response approximation. In both papers [8l[2Tj. the coupling strength among the 
elements of the chain is suggested as a suitable parameter besides the noise strength to efficiently 
perform desired operations. In this work, rather than dealing with a globally coupled network, 
we will consider as in [8] that the bistable units have nearest neighbor interactions. We will 
still be interested in finite arrays. While in [8] the authors focus on the response of a single 
oscillator, we will concentrate here as in 9, 2 1] on the behavior of a collective or global variable 
characterizing the entire array rather than on a single individual. By contrast with the range 
of parameters analyzed in |21j , we will be dealing here with situations where a linear response 
approximation is not valid. 

The rest of the paper is as follows. In Section[2j we introduce the model system, fix notation 
and indicate the relevant variables that we will use to quantify the SR effect. In Section [3] we 
present the main results of our numerical simulations. In Section [4] we comment on the main 
conclusions of our paper. 

2 The model system 

We consider a set of N identical bistable elements characterized by the variables Xi(t) (i = 
l,...,N) with nearest neighbor interactions. The dynamics is given by stochastic evolution 
equations (in dimensionless form) of the type 



subject to the conditions XN+i(t) = Xi(t), xa(t) — xpf(t). The external driving force is periodic 
in time with period T, i. e., F(t) = F(t + T). The term £,(t) represents a white noise with zero 
average and (£i(^)£j ( s )) = 6ijS(t— s). In the N — > oo limit, this model becomes the classical 4> A 
model analyzed in [22]. In the absence of driving, it has also been used to model the dynamics 
of voltage pulses along myelinated nerves [55]. In the context of the ratchet effect, a similar 
model has been studied in [24] , 

We define a collective variable S(t) as 



We will concentrate on the SR effects associated with the collective variable when the system 
size, N, is kept finite and the amplitude of the driving term is small. By small we mean that 
when the driving force acts on a single isolated unit, SR is well described by linear response 
theory. 

We will use the signal-to-noise ratio (SNR) of the collective behavior as the quantifier of the 
SR effects. Its definition requires the evaluation of the one-time correlation function defined as 



±i{t) = Xi(t) - xf{t) + -[xi- x {t) + x i+1 (t) - 2xi(t)} + V2D&(t) + F(t), 



(1) 




(2) 



3=1 




(3) 
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The notation (. . .) indicates an average over the noise realizations and the subindex oo indicates 
the long time limit of the noise average, i. e., its value after waiting for t long enough for the 
transients to die out. As indicated in our previous work [9], we can write 

L(t) = L C oh(r) + L incoh (T), (4) 

where the coherent part, L con (T), 



1 



T 



icoh(r) = - / dt <5(t)) 00 (5(t + t))oo, (5) 
1 Jo 

is periodic in t with the period of the driving force, while the incoherent part, Lincohtr) arising 
from the fluctuations of the output S(t) around its average value, decays to zero as r increases. 
The output SNR, R out , is 

Rout = lim — = — — = -= — , (6) 

where Q = 2n/T is the fundamental frequency of the driving force Fit), L con (f2) is the cor- 
responding Fourier coefficient in the Fourier series expansion of i C oh(T), and £; nco h(J?) is the 
Fourier transform at frequency Q of L- lnco h(T). 

We will also discuss the SR gain, G, defined as [9] 

G = §^, (7) 

where i?i np is the SNR of the random input process formed by the arithmetic mean of the 
individual noise terms plus the deterministic driving force F(t), namely, F(t) + with 
= N~ x J2i=i The gain can be seen as a dimensionless parameter that compares the 
output SNR to that of the input and, in this sense, it measures the quality of the output relative 
to the input. 



3 Results 

Following the numerical procedure detailed in our previous works [151120) . we have estimated 
the coherent and incoherent part of the collective correlation function L(t), by integrating 
the Langevin equations, Eq. (JT|) and averaging over several thousand noise realizations. With 
this information, we evaluate numerically the integrals defining the Fourier coefficients at the 
driving frequency and, using Eqs. © and ([7]), the collective SNR and gain for a wide range of 
parameter values. We have used a weak driving rectangular force given by 

F(t) = (8) 

where n(t) — [2t/T\, [z\ is the floor function of z, i.e., the greatest integer less than or equal 
to z. We have also considered the case of a sinusoidal driving term Asm(Slt). In both cases, 
the force amplitude will be taken to be A = 0.1, much smaller than the barrier height of an 
isolated bistable unit while the fundamental frequency will be Q = 0.01. 

As noticed in [5], in the case of noninteracting units (0 — 0), the SNR of the collective 
output is N times larger than that of a isolated unit driven by the same force. Nonetheless, as 
discussed in [9], the gain associated with the collective output is just the same as the one of a 
single, isolated, unit. Thus, for the weak forces that we are considering here, the collective gain 
does not exceed unity, in agreement with the predictions of the linear response theory |15 p 25 [ 
126] . The introduction of interactions between the bistable units drastically changes this picture 
and the enhancement of the SNR leads to the possibility of observing gains larger than unity 
for the weak driving forces considered. These facts can be observed in Fig. Q] 
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Fig. 1. The collective SNR, i? ut, (upper panel) and the collective gain, G, (lower panel) vs. the noise 
strength D for an array of TV = 10 bistable units driven by a rectangular driving force with amplitude 
A — 0.1, fundamental frequency Q = 0.01 and several values of the coupling parameter: 8 = (crosses), 
0.2 (circles), 0.5 (diamonds), 1.0 (squares), and 1.5 (triangles). Lines are a guide to the eye. 




t 



Fig. 2. The coherent part, L co h{t), (left panel) and the incoherent part, Li nco h(t), (right panel) of 
the correlation function of the collective variable for several values of the coupling parameter: 9 — 
(solid line), 0.5 (dotted), 1.0 (dashed), and 1.5 (dot-dashed), corresponding to the noise strength values 
D = 0.1, 0.16, 0.24, and 0.28, respectively. These noise values correspond to the peaks observed in i? ou t 
in the upper panel of Fig. [1] Other parameter values: TV = 10, amplitude A = 0.1 and fundamental 
frequency SI = 0.01. 



The non-monotonic behavior of the SNR for the collective variable with the noise strength 
observed in the upper panel of Fig. [T] is indicative of the SR phenomenon. In the upper panel of 
Fig. [TJ we depict the behavior of the global i? ut for an array of TV — 10 identical particles with 
nearest neighbor coupling. The values of the interaction parameter range from small values 
(6 = 0.2) to rather large ones (6 = 1.5). The i? ut peak value depends on the coupling strength 
in such a way that as 9 is increased, the noise value at which i? ou t reaches its maximum is 
shifted slightly to higher values. 

It is interesting to compare the time behavior of the coherent, L co h(t), and incoherent, 
Lincoh(t), parts of the collective correlation function. In the left panel of Fig. [21 we depict the 
time behavior of the coherent part for several values of the interaction strength 0. For each value 
of 6, the noise strength value is that at which the SNR is maximal. The periodicity of L C oh(£) 
is clearly demonstrated. Its amplitude is only slightly dependent on the interaction strength. 
Its Fourier component at the fundamental driving frequency is precisely the numerator of i? ut- 
Their peak values are not much different from those obtained in a single bistable unit driven 
by the same force at the same noise strength. 

In the right panel of Fig. [21 the behaviors of the corresponding incoherent parts are de- 
picted. It is remarkable the fast decay of the fluctuations as well as their small values. As the 
denominator of the SNR is the Fourier component of those decaying functions at the driving 
frequency, it is clear that those contributions are small. Consequently, the values of the SNR are 
expected to be much enhanced with respect to those values typical of SR in the linear regime. 
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Fig. 3. Solid line: The incoherent part of the correlation function of the global variable for an array of 
N = 10 bistable units driven by a rectangular driving force with amplitude A = 0.1 and fundamental 
frequency S7 = 0.01. Dashed line: the equilibrium time correlation function of the global variable for 
the same array in the absence of driving. In both cases the coupling parameter is 9 — 0.5 and the noise 
strength value D — 0.16. 

The enhancement of SR effects in arrays of interacting bistable units with respect to those in 
individual units is then basically a consequence of the strong reduction of the fluctuation level 
with respect to the one found in an single driven unit. 

It should also be noted in the lower panel of Fig. [l] that the gain can be larger than unity for 
some ranges of noise strength values. This feature is a clear indication that the SR phenomenon 
observed in the array for the parameter values considered can not be described within the limits 
of a linear response theory. To further understand why linear response theory fails in the cases 
considered here, it seems useful to compare the behavior of the incoherent part of the one-time 
correlation function of the global variable and that of the equilibrium correlation function of the 
same global variable in an un-driven system. An example of such comparison is depicted in Fig. 
[3] for a coupling strength 9 = 0.5. The graph clearly indicates that the equilibrium fluctuations 
are much larger and longer lasting than the fluctuations about the average behavior in the driven 
system. In the linear response theory description of SR, (see, for instance, [1] and references 
therein), it is assumed that the correlation function of the fluctuations around the average 
behavior in a driven system can be safely approximated by their corresponding equilibrium 
values in a un-driven one. This is clearly not the case for the the system at hand. A detailed 
study of the validity conditions of linear response theory can be found in [7flT25] . 

The peak values of i? ut also show a non-monotonic behavior with 9 as depicted in Fig. |4] 
As 9 is raised from zero up to 0.5, there is an increase in the peak of the SNR value. This is due 
to the combination of two effects: i) the faster decay of the correlation function as 9 is increased 
(see right panel in Fig. [2|) with the subsequent decrease of the denominator in the SNR, and ii) 
the larger amplitude of the coherent part (see left panel in Fig. [2]). On the other hand, as the 
9 values is further increased, the contribution of the coherent part decreases, while that of the 
incoherent part increases and, consequently, the SNR decreases as the coupling term increases. 
The gain also shows a non-monotonic behavior with 9 as depicted in Fig. [4] 

The complexity of the TV-dimensional potential surface where the interacting particles move 
renders the task to give a simple explanation of the observed non-monotonic behavior with 9 
a difficult one. For an A-dimensional surface in the absence of interactions and driving forces, 
the whole surface is symmetrical about the origin with barriers of equal heights along each axis. 
The observed non-monotonic behavior can be rationalized in terms of the periodic rocking of 
the bistable potential independently along each axis. Each particle jumps over its corresponding 
barrier independently of the other particles, under the influence of the driving force and the 
noise. The central limit theorem can be safely used for independent subunits, so that the increase 
in the SNR values is just a size effect. 

On the other hand, for coupled systems, the individual stochastic processes Xi(t) are no 
longer independent and the central limit theorem alone is not enough to understand the reported 
results. The deformation of the potential surface due to the external driving and the interaction 
term might very well reduce the height of the barriers, eliminate some of them and alter the 
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Fig. 4. The peak values of R ou t and G vs. the coupling strength 6 for an array of N = 10 elements 
driven by a rectangular force with A — 0.1 and fundamental frequency Q — 0.01. 



location of the minima. Then, one can not rule out the possibility of the existence of new 
paths facilitating the transitions between the attractors. This being the case, a reduction of the 
fluctuation levels besides the one coming from the system size is to be expected. For a fixed 
small driving amplitude, the amount of distortion of the potential relief must depend on the 
strength of the coupling term 9. For very small values of 9, particles located beyond nearest 
neighbor positions are expected to be weakly correlated, with the correlation length increasing 
as 9 increases. One might expect an increase on the SNR values as 9 increases from zero, as 
the incoherent part of the correlation function basically decreases. This is due to the increasing 
easiness of transitions between attractors along the new paths. As 9 is further increased, the 
distortion of the potential relief will increase, but at the same time, most of the particles along 
the chain will start to be strongly correlated. In other words, the motion of the collective variable 
will resemble more and more the motion of a particle on a single bistable potential. The chain 
will behave more and more like a rigid object. Jumps over the barriers become increasingly 
more difficult and this leads to a decrease of the SNR values. There must be an intermediate 
interaction strength value so that a maximum SNR of the output is achieved. 

It is interesting to study what happens in the case of a sinusoidal driving force with the 
same amplitude A = 0.1 and frequency fl = 0.01. In Fig. [5] we depict R out (lower panel) and 
G (upper panel) vs. D for several values of the coupling strength. As it can be seen, the SNR 
values are substantially smaller than the ones observed for a rectangular driving force (compare 
with Fig. [I]). The gain is always below unity by contrast with the rectangular driving force, 
where the gain can reach values larger than 1. 

In Fig. [6]we show the time dependence of the coherent and incoherent parts of the correlation 
function for an array of N — 10 particles driven by the sinusoidal force. By comparison with 
Fig. [2] we see that the reason why the SNR values are much smaller for a sinusoidal driving 
than for a rectangular one is mainly that £mcoh(i) is larger in the single frequency case than 
in the multi-frequency one. It is not correct to conclude that because the gain in the case of 
sinusoidal driving is less than unity, a linear response theory description is adequate. As shown 
in Fig. [Jj Lincoh^) for a sinusoidal driving is very different from the equilibrium correlation 
function in an un-driven system. 

A non-monotonic behavior of the SNR with the coupling constant 9 for the sinusoidal driving 
also exists. The qualitative explanation given before for the rectangular driving still stands. 
Note, nonetheless, that the rocking of the multi-dimensional energy surface brought up by the 
sinusoidal driving is less drastic than the one produced by the rectangular one. Even though 
the forces have the same amplitude and fundamental frequency, the sinusoidal force introduces 
a bias in the energy relief with respect to that in the zero force case which is continuously 
changing with time. On the other hand, the rectangular signal keeps the surface biased most of 
the time, except during its instantaneous changes of value. The system has ample time to relax 
during those biased intervals and this explains why the fluctuations are so drastically reduced 
in the case of a rectangular driving. 
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Fig. 6. The time behavior of L C oh and Z/incoh for several values of the coupling strength for an array 

= 0.1 and frequency Q — 0.01. 
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Fig. 7. Comparison of the time behavior of Lincoh in an array of N = 10 at equilibrium or driven by 
either a sinusoidal force with A — 0.1 and frequency £2 = 0.01, or a rectangular one with the same 
amplitude and fundamental frequency. Other parameter values: 9 — 0.5; D = 0.16 



4 Conclusion 

In this work, we have explored the phenomenon of SR in finite arrays of noisy bistable units 
with nearest neighbor coupling, driven by time periodic forces of weak amplitude. Rather than 
analyzing the modification of a single unit behavior with respect to the one in the absence of 
any interaction, we have focused our attention in a collective variable describing the dynamics 
of the array as a whole. 

A strong enhancement in the SR effects associated with the collective variable with respect to 
the one found in a system formed by a single unit has been demonstrated by means of numerical 
simulations. In particular, we have shown that, even for rather weak input amplitudes, the SNR 
has a non-monotonic behavior with the noise strength. The different values of the SNR quantifier 
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are much larger than the corresponding ones found in a single unit system. Furthermore, the 
SR gain reaches values higher than unity for some values of the parameters, indicating that the 
array is operating in a nonlinear regime well beyond the limits of the regimes described by linear 
response theory. This is to some extent surprising as, for the driving amplitudes and frequencies 
considered, the response of a single unit system is well described by the linear response theory. 

There are in principle two main reasons for the enhanced effects reported here. On the one 
hand, based on the central limit theorem, one can expect a decrease in the fluctuation levels 
with respect to those found in single unit systems simply because of the system size. On the 
other hand, as we have demonstrated in our numerical simulations, the size effect mechanism 
is not enough to explain the reported behavior. The presence of coupling terms between the 
subsystems is an essential ingredient to obtain gains larger than unity. 

We believe that the main reason for the strong enhancement observed is due to the drastic 
reduction of the value of the incoherent correlation function and on its correlation time induced 
by the system size, the external driving and the coupling term between the different sub-units. 
Indeed, the comparison of the incoherent correlation function with the equilibrium correlation 
function in finite systems with the same size and coupling parameter, shows the relevance of 
the driving force. 

An enhancement of the SR effects in finite size arrays with global coupling (mean field 
coupling) have also been reported by us. One of the goals of the present work is to demonstrate 
the robustness of our results, regardless of the topology of the connections between the different 
subunits. Indeed, a comparison of the numerical results reported here with those in [20] indicate 
that, at least for small size arrays (N = 10), the i? ou t and gain values are not much different. 

We have also noted that there exists a non-monotonic behavior of the SNR with respect to 
the coupling parameter 9, besides the usual non-monotonic behavior with the noise strength. 
There are some indications that non-monotonic behaviors with respect to the system size might 
also exist. Indeed, system size resonances have been analyzed by other groups for systems of 
globally coupled nonlinear oscillators using cumulant expansion techniques |29j or nonequi- 
librium potentials [30j . Within a linear response theory description, they report system size 
resonance effects. We are presently investigating the issue of the dependence on the system 
size in coupled arrays (global and local coupling) for parameter regimes well beyond the linear 
response theory limits. 
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